On the Energy Spectra of GeV/TeV Cosmic Ray Leptons

Recent observations of cosmic ray electrons from several instruments have revealed various degrees of deviation in the measured electron energy distribution from a simple power-law, in a form of an excess around 0.1 to 1 TeV energies. An even more prominent deviation and excess has been observed in the fraction of cosmic ray positrons around 10 and 100 GeV energies. These observations have received considerable attention and many theoretical models have been proposed to explain them. The models rely on either dark matter annihilation/decay or specific nearby astrophysical sources, and involve several additional assumptions regarding the dark matter distribution or particle acceleration. In this paper we show that the observed excesses in the electron spectrum may be easily reproduced without invoking any unusual sources other than the general diffuse Galactic components of cosmic rays. The model presented here assumes a power-law injection of electrons (and protons) by supernova remnants, and evaluates their expected energy spectrum based on a simple kinetic equation describing the propagation of charged particles in the interstellar medium. The primary physical effect involved is the Klein-Nishina suppression of the electron cooling rate around TeV energies. With a very reasonable choice of the model parameters characterizing the local interstellar medium, we can reproduce the most recent observations by Fermi and HESS experiments. Interestingly, in our model the injection spectral index of cosmic ray electrons becomes comparable to, or even equal to that of cosmic ray protons. The Klein-Nishina effect may also affect the propagation of the secondary e{sup {+-}} pairs, and therefore modify the cosmic ray positron-to-electron ratio. We have explored this possibility by considering two mechanisms for production of e{sup {+-}} pairs within the Galaxy. The first is due to the decay of {pi}{sup {+-}}'s produced by interaction of cosmic ray nuclei with ambient protons. The second source discussed here is due to the annihilation of the diffuse Galactic {gamma}-rays on the stellar photon field. We find that high positron fraction increasing with energy, as claimed by the PAMELA experiment, cannot be explained in our model with the conservative set of the model parameters. We are able, however, to reproduce the PAMELA (as well as Fermi and HESS) results assuming high values of the starlight and interstellar gas densities, which would be more appropriate for vicinities of supernova remnants. A possible solution to this problem may be that cosmic rays undergo most of their interactions near their sources due to the efficient trapping in the far upstream of supernova shocks by self-generated, cosmic ray-driven turbulence

Upper Bound on the First Star Formation History

Our understanding of the nature of the extragalactic background light (EBL) has improved with the recent development of gamma-ray observation techniques. An open subject in the context of the EBL is the reionization epoch, which is an important probe of the formation history of first stars, the so-called Population III (Pop III) stars. Although the mechanisms for the formation of Pop III stars are rather well understood on theoretical grounds, their formation history is still veiled in mystery because of their faintness. To shed light into this matter, we study jointly the gamma-ray opacity of distant objects and the reionization constraints from studies of intergalactic gas. By combining these studies, we obtain a sensitive upper bound on the Pop III star formation rate density as $\dot\rho_{*}(z)<0.01[(1+z)/{(1+7.0)}]^{3.4}({f_{\rm esc}}/{0.2})^{-1}({C}/{3.0})\ {\rm M}_{\odot} {\rm yr}^{-1}\ {\rm Mpc}^{-3}$ at $z\ge7$, where $f_{\rm esc}$ and $C$ are the escape fraction of ionizing photons from galaxies and the clumping factor of the intergalactic hydrogen gas. This limit is a $\sim10$ times tighter constraint compared with previous studies that take into account gamma-ray opacity constraints only. Even if we do not include the current gamma-ray constraints, the results do not change. This is because the detected gamma-ray sources are still at $z\le4.35$ where the reionization has already finished.Comment: 5 pages, 2 figures, accepted for publication in ApJ

On the Radio and Optical Luminosity Evolution of Quasars

We calculate simultaneously the radio and optical luminosity evolutions of quasars, and the distribution in radio loudness R defined as the ratio of radio and optical luminosities, using a flux limited data set containing 636 quasars with radio and optical fluxes from White et al. We first note that when dealing with multivariate data it is imperative to first determine the true correlations among the variables, not those introduced by the observational selection effects, before obtaining the individual distributions of the variables. We use the methods developed by Efron and Petrosian which are designed to obtain unbiased correlations, distributions, and evolution with redshift from a data set truncated due to observational biases. It is found that the population of quasars exhibits strong positive correlation between the radio and optical luminosities. With this correlation, whether intrinsic or observationally induced accounted for, we find that there is a strong luminosity evolution with redshift in both wavebands, with significantly higher radio than optical evolution. We also construct the local radio and optical luminosity functions and the density evolution. Finally, we consider the distribution of the radio loudness parameter R obtained from careful treatment of the selection effects and luminosity evolutions with that obtained from the raw data without such considerations. We find a significant difference between the two distributions and no clear sign of bi-modality in the true distribution for the range of R values considered. Our results indicate therefore, somewhat surprisingly, that there is no critical switch in the efficiency of the production of disk outflows/jets between very radio quiet and very radio loud quasars, but rather a smooth transition. Also, this efficiency seems higher for the high-redshift and more luminous sources in the considered sample.Comment: 15 pages, 15 figures, accepted to ApJ, updated to in press versio

A Multi-Code Analysis Toolkit for Astrophysical Simulation Data

The analysis of complex multiphysics astrophysical simulations presents a unique and rapidly growing set of challenges: reproducibility, parallelization, and vast increases in data size and complexity chief among them. In order to meet these challenges, and in order to open up new avenues for collaboration between users of multiple simulation platforms, we present yt (available at http://yt.enzotools.org/), an open source, community-developed astrophysical analysis and visualization toolkit. Analysis and visualization with yt are oriented around physically relevant quantities rather than quantities native to astrophysical simulation codes. While originally designed for handling Enzo's structure adaptive mesh refinement (AMR) data, yt has been extended to work with several different simulation methods and simulation codes including Orion, RAMSES, and FLASH. We report on its methods for reading, handling, and visualizing data, including projections, multivariate volume rendering, multi-dimensional histograms, halo finding, light cone generation and topologically-connected isocontour identification. Furthermore, we discuss the underlying algorithms yt uses for processing and visualizing data, and its mechanisms for parallelization of analysis tasks.Comment: 18 pages, 6 figures, emulateapj format. Resubmitted to Astrophysical Journal Supplement Series with revisions from referee. yt can be found at http://yt.enzotools.org

The VLBA Imaging and Polarimetry Survey at 5 GHz

We present the first results of the VLBA Imaging and Polarimetry Survey (VIPS), a 5 GHz VLBI survey of 1,127 sources with flat radio spectra. Through automated data reduction and imaging routines, we have produced publicly available I, Q, and U images and have detected polarized flux density from 37% of the sources. We have also developed an algorithm to use each source's I image to automatically classify it as a point-like source, a core-jet, a compact symmetric object (CSO) candidate, or a complex source. The mean ratio of the polarized to total 5 GHz flux density for VIPS sources with detected polarized flux density ranges from 1% to 20% with a median value of about 5%. We have also found significant evidence that the directions of the jets in core-jet systems tend to be perpendicular to the electric vector position angles (EVPAs). The data is consistent with a scenario in which ~24% of the polarized core-jets have EVPAs that are anti-aligned with the directions of their jet components and which have a substantial amount of Faraday rotation. In addition to these initial results, plans for future follow-up observations are discussed.Comment: 36 pages, 3 tables, 13 figures; accepted for publication in Ap

Flux and Photon Spectral Index Distributions of Fermi-LAT Blazars and Contribution to the Extragalactic Gamma-ray Background

We present a determination of the distributions of gamma-ray flux - the so called LogN-LogS relation - and photon spectral index for the 352 blazars detected with a greater than approximately seven sigma detection threshold and located above {+-} 20{sup o} Galactic latitude by the Large Area Telescope of the Fermi Gamma-ray Space Telescope in its first year catalog. Because the flux detection threshold depends on the photon index, the observed raw distributions do not provide the true LogN-LogS counts or the true distribution of the photon index. We use the non-parametric methods developed by Efron and Petrosian to reconstruct the intrinsic distributions from the observed ones which account for the data truncations introduced by observational bias and includes the effects of the possible correlation among the two variables. We demonstrate the robustness of our procedures using a simulated data set of blazars and then apply these to the real data and find that for the population as a whole the intrinsic flux distribution can be represented by a broken power law of slopes -2.37 {+-} 0.13 and -1.70 {+-} 0.26, and the intrinsic photon index distribution can be represented by a Gaussian with mean 2.41 {+-} 0.13 and 1{sigma} width of 0.25 {+-} 0.03. We also find the intrinsic distributions for the sub-populations of BL Lac and FSRQs type blazars separately. We then calculate the contribution of blazars to the diffuse cosmic gamma-ray background radiation to be 28% {+-} 19%

Generalized Attractor Points in Gauged Supergravity

The attractor mechanism governs the near-horizon geometry of extremal black holes in ungauged 4D N=2 supergravity theories and in Calabi-Yau compactifications of string theory. In this paper, we study a natural generalization of this mechanism to solutions of arbitrary 4D N=2 gauged supergravities. We define generalized attractor points as solutions of an ansatz which reduces the Einstein, gauge field, and scalar equations of motion to algebraic equations. The simplest generalized attractor geometries are characterized by non-vanishing constant anholonomy coefficients in an orthonormal frame. Basic examples include Lifshitz and Schroedinger solutions, as well as AdS and dS vacua. There is a generalized attractor potential whose critical points are the attractor points, and its extremization explains the algebraic nature of the equations governing both supersymmetric and non-supersymmetric attractors